% Devin Koepl

function [y, kx, ky, bx, by] = dissipative_environment
y   = @(x) zeros(size(x));


kx  = @(x) zeros(size(x));
ky  = @(x) zeros(size(x));

%     uniform_damping = 1e3;

uniform_damping = evalin('base','uniform_damping');

bx  = @(x) uniform_damping*ones(size(x));
by  = @(x) uniform_damping*ones(size(x));
end
